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^ ' Aging phenomena of short-range Ising spin glass models have been investigated using 

M . Monte Carlo simulations. It is found that in the low-temperature spin-glass phase the mean 

L ■ domain size exhibits a crossover from a power-law growth associated with the critical fluctu- 

• i—j ' ation at the transition temperature to slower growth inherent in the low-temperature phase. 

' O \ The temperature dependence of the growth law of the domain size can be almost explained by 

■^—i , this crossover. We also find that the spin-autocorrelation function in the quasi-equilibrium 

2 ■ regime follows the expected scaling behavior from the droplet theory expressed in terms of 

S ' the mean domain size and a characteristic length scale of droplet excitations. 

^'. 

Ch ' §!• Introduction 

o\ 

I— 1 1 Spin glasses (SG) exhibit characteristic slow dynamics below the SG transition 

temperature. ^>^> Recently such slow dynamics, in particular aging phenomena, has 

^ , been much attractive both in experimental and theoretical studies. Several attempts 

■^ I have been made so far to explain the slow dynamics. There are mainly two distinct 

■ ■ phenomenological approaches; one is a phase-space approach, in which the dynamics 

^^ , is described by a diffusion in the hierarchically constructed phase space, inspired 

^~~i \ by the mean-field theory which suggests multi-valley structure of the phase space. 

^Z^ • The other is a real-space picture based on the scaling theory^ in which low- lying 

excitations are attributed to connected clusters reversed from one of two ground 

C^ ' states. There has been, however, no satisfactory description from a microscopic 

^ ■ model, except for the dynamical mean-field theory. ^' In this paper we present the 

results on the non-equilibrium dynamics obtained by large-scale Monte Carlo (MC) 

^ ' simulations on short-range Edwards- Anderson (EA) Ising SG models. Our analyses 

O ■ of the obtained data are based on the droplet theory. We believe that, even when the 

phase-space approach may give us a correct description for the slow dynamics, what 

.^ ' really occurs in the real space is also indispensable for thorough understandings of 

^ ■ the aging phenomena. 
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H I Let us explain briefly the droplet theory. ^^ According to the theory, aging process 

is described by coarsening of domain walls, which is driven by successive flipping of 
thermally activated droplets. During isothermal aging up to waiting time tw after 
quench, domains with the mean size -R(t«,) separating different pure states have 
grown up. Within each domain, small droplets of size L ~ L{t) <^ R{tw) are 
thermally fluctuating within a time scale of r as in equilibrium. The typical value 
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of their excitation gap F^"^ scales as 

Ff P ~ r(L/Lo)^ (M) 

where T is the stiffness constant and Lq is a microscopic length scale, and that of 
free-energy barrier B^^ also scales as 

B'l^ ~ A{L/L,)\ (1-2) 

where Z\ is a characteristic free-energy scale. 

As compared with the equilibrium, some droplets which touch the domain wall 
with the length scale R{tw) could reduce their excitation gap from (1-1). This effect 
is estimated by the droplet theory as the reduction of the averaged excitation gap 
which is given by 

F^P=reff(L/Lor, (1-3) 

with the effective stiffness constant 

%^ = r(l-c,{L/Rf-'), (1-4) 



where c„ is a numerical constant. Physical quantities such as the spin autocorrelation 
function are estimated in terms of such length scales by taking into account statistical 
weights of the droplet excitations appropriately. The growth law of the length scales 
of the domain R(tw) and the droplet L{t) is given by 

i?(t),L(t)~(^-ln(t/To)j , (1-5) 

where tq is a microscopic time scale. One notices that the droplet theory for aging 
phenomena consists of two almost independent steps; the scaling argument based 
on the typical length scales of i?(t«,) and L{t) and the growth law of these length 
scales. The main purpose of the present work is to test these two steps separately. 

The present paper is organized as follows: in the next section after introducing 
the model system studied, time evolution of the length scale of domain wall is dis- 
cussed. The results of spin- autocorrelation function in the quasi-equilibrium regime 
of isothermal aging are presented in Sect. 4. 

§2. Model and Method 

We focus on the four-dimensional (4D) Ising SG model, because its static critical 
properties have been established in the sense that a SG phase transition occurs at 
finite temperature with a rigid order parameter. The 4D EA Ising SG model is 
defined by Hamiltonian, 

n = -J2J^,S^Sj, (2-1) 

where the sum runs over nearest-neighbor sites and the Ising variables Si are defined 
on a hypercubic lattice with periodic boundary conditions. The interactions are 
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bimodal variables (Jij = ±1) distributed randomly with equal probability. The 
simulation method is the standard single-spin-flip Monte Carlo (MC) method using 
two-sublattice dynamics with the heat-bath transition probability. Using the MC 
method, we have simulated aging processes after a rapid quench from T = oo to the 
SG phase. The system size studied is mainly L = 24, while the size L = 32 is partly 
studied. We have found no significant difference in the data of L = 24 and 32. 

§3. Growth Law of Domain Size 

In order to extract a length scale characterizing the growth of ordering, we 
calculate spatial replica correlation function in off-equilibrium, defined as 

G{r,t) = J2ist\t)sf\t)slll{t)sfJ^{t)), (3-1) 

i 

where a and (3 denote the replica indices which are updated independently and 
with different initial spin configurations. The bracket (• • •) denotes the average 
over independent bond realizations. In our simulations, only one MC sequence is 
performed for each random bond configuration. We extract the mean domain size 
R{t) by directly fitting G{r, t) to an exponential form. 

In Fig. 1 we show time dependence of R{t) at the SG transition temperature 
Tc(~ 2.0J) and below. Just at Tc, the length scale is expected to grow as a power 
law with the dynamical critical exponent z, R{t) ~ t^'^, irrespective of the physical 
picture underlying the ordered phase. As expected, it is found that the length scale 
follows a power law. The exponent z is estimated to be 4.98(5), which is consistent 
with that of the previous work. ^^ Below Tc, R{t) grows with time slower and slower 
as temperature decreases. We try to see crossover between the critical fluctuation 
and slow dynamics inherent in the low-temperature phase. In off-equilibrium both at 
the critical and the off-critical temperatures, the length scale would exhibit a power 
law in short length and time regime where the critical fluctuation dominates the 
dynamics. We assume that such microscopic length Rq and time tq are the correla- 
tion length and time associated with critical fluctuation in equilibrium, respectively. 
Thus, we propose a scaling form 

R{t)/Ro = git/To), (3-2) 

with Rq = \T — Tc\''^ and tq = \T — Tc\~'^'^. The scaling plot of R{t) is shown in 
Fig. 2. As expected from a standard scaling theory of critical phenomena, the scaling 
function g{x) for smaller x exhibits a power law x^'^ associated with the critical 
temperature. We find a significant deviation from the power law at longer times, 
suggesting that the characteristic slow dynamics of the SG phase takes place there. 
In fact, the functional form is not incompatible with a power law of ln(t) predicted 
by the droplet theory (1-5). It is noted that the strong temperature dependence of 
R{t) shown in Fig. 1 can be almost explained by introducing the microscopic units 
Ro and tq associated with the critical fluctuation, while it is not sure whether the 
asymptotic form at longer times is also described by a universal scaling function or 
not. 



K. Hukushima, H. Yoshino and H. Takayama 



0.5 
0.4 





.■' 


z~4.98(5) 


j?i/' ' 








:,.^ ° 




■ ■ 1 


.■ii'&. i • 


r//=2.o ^-m-^ 
1.8 
1.6 


' i 


1.2 - 




0.8 



1000 10000 100000 

t 



Fig. 1. R{t) of the 4D Ising SG model as a 
function of f at Tc and below Tc. 
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Fig. 2. Scaling plot of i?(f ) of the 4D Ising SG 
model, where the SG transition tempera- 
ture Tc = 2.0J and the critical exponent 
ly = 0.93 are fixed. 



§4. Scaling analysis in quasi-equilibrium 

Recent studies have revealed that the off-equihbrium dynamics in the SG is sepa- 
rated into two characteristic time regimes. One is a short-time regime, called "quasi- 
equilibrium regime" , and the other is a long-time regime, called "aging regime" . A 
typical observable to see such two time regimes is the spin auto-correlation function 



C(T;t, 



-E 



{Si{tyj)Si{T + t^)), 



(4-1) 



where t^ denotes a waiting time after the rapid quench. Our interest is in its behavior 
in the quasi-equilibrium regime, namely t„, ^ r. Based on the droplet argument 
using the effective stiffness constant (1-4), the behavior of C{t; tw) is explicitly given 

by 7) 



C(r;tw)=Ceq(r) + 



fL{T) 



with the equilibrium part Ccq(r)^^ 

A 



d-e 



+ 



Ceq(T) = g'EA 



{L{t)/L,)0' 



(4-2) 



(4-3) 



where ^ea is the EA order parameter and c and A are numerical constants. This 
gives us an extrapolation form of the large r limit, namely a way of determining the 
EA order parameter. An empirical form C{t) = qea + a/t'^ has been used frequently 
for estimating qea from the spin autocorrelation function. ^' This is true only if the 
time dependence of length scale L{t) in (4-3) is a power law. As seen in the last 
section, however, the observed length scale R{t) in this model exhibits the crossover 
from the critical power law to slower growth at large t. 

We observe the autocorrelation function C{T;tw) at T/J = 1.2 well below Tc 
and check the scaling form (4-2) by making use of the R{t) estimated through G(r, t) 
in the last section. According to the droplet theory, both R{tyj) and L{t) exhibit 
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the same time dependence. Also, the microscopic length scale Lq is assumed to be 
the same as Rq. For fixed r, the autocorrelation function is expressed as a function 
of R{tw). Using the estimated value of 6{= 0.82), ^^ a simple linear fitting gives 
the equilibrium autocorrelation function Ceq(T) in the large t^ limit. Collecting 
Ceqij) for each r thus extracted, we confirm directly the droplet prediction (4-3) as 
shown in Fig 3 and determine the equilibrium EA order parameter. The value of 
gEA estimated to be 0.58(1) is compatible with the recent estimation from static MC 
simulation. ^^' 

Next we discuss correction to the equilibrium limit. The expression (4-2) suggests 
that the correction term AC{t; t^) = C{t; t^) — Ccq(r) multiplied by L^{t) becomes 
only a function of L{T)/R{tyj). As shown in Fig. 4, we confirm this scaling predic- 
tion. For the limit of L{t) / R{tw) <^ 1, the scaling function shows (L(r)/i?(t^)) ~ 
consistent with (4-2). 
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Fig. 3. Equilibrium autocorrelation function 
extracted to the large t^ limit. The line 
represents a fitting according to the expres- 
sion (4-3). 
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Fig. 4. Scaling plot of the correction term of 
C{r;t^). The line has an expected slope 
d — 6 for small L{t) / R{tw) limit. 



§5. Discussion and Summary 



Let us compare the present results with those obtained in 3D Ising SG models. 
In three dimensions, results on the growth law by numerical simulations ^^^ " ^^-^ as 
well as experiments ^'^> are well fitted to a power law as 



R{t) ~ t 



l/z{T) 



(5-1) 



where the exponent l/z{T) is proportional to temperature T and continuously con- 
nects with the dynamical critical exponent zatTc- It is not clear yet how to interpret 
physically such a power law with temperature-dependent exponent l/z{T). One of 
the possibilities is the crossover observed in the preset work on 4D Ising SG model. 
Certainly it is worth examining the crossover effect of the critical fiuctuation in the 
3D model. 

On the other hand, the scaling argument in terms of the length scales i?(t«,) 
and L(t) successfully explains the aging behavior of the correlation function in the 
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quasi-equilibrium regime also in three dimensions. Recently it has been confirmed 
not only in the isothermal but also in temperature-shift aging processes, ^^-* which 
are basic experimental procedures frequently used. 

To conclude, we have investigated non-equilibrium dynamics after the temper- 
ature quench from infinity to the SG phase in the 4D Ising SG model using Monte 
Carlo simulations. We have studied the growth law of the mean domain size by ana- 
lyzing time evolution of the spatial replica correlation functions. We have found that 
the growth law shows a crossover from the critical regime to the low-temperature 
one and its main temperature dependence within time range of our simulation can 
be explained by this crossover. 

We have also analyzed the spin autocorrelation function in the quasi-equilibrium 
regime. The off-equilibrium correction of the correlation function to its equilibrium 
limit, namely, violation of the time translational invariance, in the quasi-equilibrium 
limit can be explained by the scaling argument in terms of the characteristic length 
scales R{tw) and I/(t), as observed already in the 3D Ising SG model. ^^ These results 
strongly suggest that the analysis consisting of the two almost independent steps is 
promising for understanding the aging phenomena in low-dimensional SG systems. 
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